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ABSTRACT 

Orbits in triaxial ellipsoidal potentials precess about either the major or minor axis of 
the ellipsoid. In standard perturbation theory it can be shown that a circular orbit will 
precess about the minor axis if its angular momentum vector lies in a region bounded 
by two great circles which pass through the intermediate axis and which are inclined 
with minimum separation it from the minor axis, where it = arctan((£? 2 — C 2 )/(A 2 — 

1/2 

B 2 )) and A, B and C are the axis ratios, A > B > C . We test the accuracy of 
this formula by performing orbit integrations to determine is, the simulated turnover 
angle corresponding to it- 

We reach two principal conclusions: (i) is is usually greater than it, by as much 
as 12 degrees even for moderate triaxialities, A/ 1.2 < B < C/0.8. This reduces the 
expected frequency of polar rings, (ii) is is not a single, well-defined number but can 
vary by a few degrees depending upon the initial phase of the orbit. This means that 
there is a reasonable probability for capture of gas onto orbits which precess about 
both axes. Interactions can then lead to substantial loss of angular momentum and 
subsequent infall to the galactic centre. 
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Polar ring galaxies are systems in which a ring of gas (and/or 
young stars) is seen to orbit about the major axis of an early- 
type galaxy. In many cases it has been established that the 
host galaxy is rotating at right-angles to the ring and this 
is sometimes taken to be part of the definition of a polar 
ring system. There are just 6 confirmed polar rings but 27 
good candidates and many more possibles — see Whitmore 
et al. (1990) for a review. 

The host galaxy in polar ring systems often appears to 
be an SO but this leads to theoretical problems: in an oblate 
axisymmetric system all orbits will precess about the minor 
axis at a rate which is a function of radius (typically the pe- 
riod is proportional to radius). This differential precession 
will cause the polar ring to fragment. This can be overcome 
if a small degree of triaxiality is assumed as orbits whose an- 
gular momentum is sufficiently close to the major axis will 
then precess around the major rather than the minor axis. 
In general accreted discs of gas will have an arbitrary orien- 
tation. In this case differential precession about a symmetry 
axis will lead to dissipation and collapse into the plane per- 
pendicular to that axis. 

The regions of phase space which lead to precession 
about the major or minor axes can be investigated using per- 



turbation theory in Hamiltonian mechanics. It is assumed 
that the unperturbed potential is spherically symmetric and 
that the orbits are circular. The Hamilton is then expanded 
in spherical harmonics and the time-averaged perturbing 
potential calculated (see, for example, Steiman-Cameron & 
Durisen, 1984). This results in an expression of the form 



< $i >oc C 20 (3 sin 2 i-2) + 6C22 cos 2D. sin 2 



(1) 



where Q and i are the node and inclination of the orbit, 
as illustrated in Figure [l], and G20 and C22 are constants. 
Orbits will precess along lines < $1 >=constant as shown in 
Figure^. If the angular momentum of the orbit, J, lies within 
the shaded region then it will precess about the major axis, 
otherwise it will precess about the minor axis. The dividing 
lines between the two regions are great circles which pass 
through the intermediate axis and which are inclined at an 
angle 



it 



. ( C20 + 2C22 \ ' 
arCSm Uo-2C722 J 



(2) 



to the minor axis. Hence the probability that a randomly 
inclined disk will precess around the major axis is 

/=1-^. 

TV 
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2 INTEGRATION METHOD AND RESULTS 



2.1 The integration method 

The numerical code used in this paper is that described by 
Pearce & Thomas (1991). It is a simple predictor-corrector 
method integrates orbits to high accuracy and we refer the 
reader to this paper for details. 
We use a potential 

(222 
x y z 
A^ + h + CP 

where A > B > C and without loss of generality we can take 
5 = 1. The time-averaged first-order perturbation potential 
takes the form of Equation ^ with 

C20 = -J7 + -gj - -^1 
C22 =\{-, 



1 



Figure 1. The orbital elements. 



which gives, on substitution in Equation 



it = arcsin 



1 

W 



1 

T 1 



arctan 



B 2 -C 2 
A 2 - B 2 



Note that this expression is different from that for a system 
with ellipsoidal density contours, p — p((x/a) 2 + (y/b) 2 + 
(z/c) 2 ), for which 



it = arcsm 



Figure 2. Precession trajectories for the angular momentum vec- 
tor of near-circular orbits in a triaxial potential. 



This formula is extensively used in the theory of polar rings. 
It is generally accepted that dissipation will cause a disk of 
gas to settle into a plane perpendicular to the axis about 
which it precesses. Even a small triaxiality can lead to a 
reasonable probability for polar ring formation. For example 
taking an ellipsoidal density distribution with axis ratios 
a = 1.01, b = 1 and c = 0.8 gives i T « 77° and / w 0.15. 

In this paper we test the validity of the perturbation 
theory by calculating orbits in an ellipsoidal triaxial poten- 
tial. We find that the measured transition angle, is, which 
divides precession about the major and minor axes, is usu- 
ally larger than the theoretical value, it- Also the boundary 
dividing the two regions of precession about the major and 
minor axes is not sharp. The precession axis depends upon 
the phase of the orbit and is can vary by a few degrees. 
The integration method and results of our simulations are 
presented in the Section 2. The conclusions are summarised 
and discussed in Section 3. 



These two formulae agree for near-spherical systems but can 
differ by 10 degrees or more for moderate triaxialities. For 
example a — 1.2, b — 1, c — 0.8 gives it — 42.1 degrees 
whereas A — 1.2, B = 1, C = 0.8 gives it = 53.6 degrees. 

Particles are given an initial tangential velocity of y/2 
which would lead to a circular orbit in a spherical poten- 
tial, A = B — C. The subsequent behaviour depends both 
upon the initial orientation of the plane of the orbit and on 
the orbital phase of the particle. These are labelled by the 
co-ordinates shown in Figure |l|. i is the inclination of the 
orbit to the x-y plane (also the angle between the angular 
momentum vector, J, and the z-axis), Q is the longitude of 
the ascending node, and x is the phase measured from the 
position of the ascending node. Because the radial velocity 
is initially zero \ corresponds also to the phase of the peri- 
or apogee of the orbit. Note, however, that this phase is not 
conserved as it would be in a Keplerian potential. 

There are too many free parameters to be able to in- 
vestigate them all in depth so we choose to restrict J to 
lie initially in the x-z plane (i.e. Q — 7r/2). If the transi- 
tion lines separating the two regions of precession about the 
major- and minor-axes are great circles (as in linear pertur- 
bation theory) then it is trivial to generalise our results to 
arbitrary J. 



2.2 Results 

The measured transition angle, is, is a function of the initial 
phase, x, as illustrated in Figure |3| for one particular choice 
of axis ratios, A — 1.1, B = 1.0 and C = 0.9. The minimum 
value of is occurs at x — and the maximum at x = t/2 
with an approximately sinusiodal variation between the two. 
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oppositely aligned. Collisions between gas clouds of similar 
mass may then reduce their velocity relative to the galac- 
tic centre to almost zero and they will be accreted into the 
core. This situation may not be as unlikely as it at first ap- 
pears because the appropriate range of is may span up to 
10 degrees even for moderate triaxialities. Also accretion of 
a gas-rich dwarf spiral or dwarf irregular galaxy is likely to 
populate a fair region of phase space and a large amount 
of dissipation is required for this to settle into a disk of gas 
clouds on circular orbits, as is often assumed for simplicity. 
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Figure 3. The measured turnover angle, %, as a function of 
initial phase for the case A = 1.1, B = 1.0 and C = 0.9. The 
theoretical value, iy, is shown as a dotted line. 

Because of this we need only present results for these two 
extremes, and the maximum and minimum transition angles 
for a range of triaxialities are given in Figure ^j. The effect 
in each case is to lower the boundary of the shaded region 
in Figure || and to smear it out over a few degrees. 

As an aside, we never see orbits which switch from pro- 
cessing about the major to the minor axis, or vice versa. In 
general the phase of the apogee of the orbit is not conserved. 
However for those orbits which lie close to the transition an- 
gle, is, the phase returns almost exactly to its original value 
after one whole precession time. We do not know why this 
is the case although it is presumably a reflection of some 
underlying conservation law. 
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3 CONCLUSIONS 

In this paper we calculate the simulated transition angle 
for precession about the major or minor axis of a triaxial 
potential and reach the following two conclusions: (i) is is 
almost always greater than it and (ii) is is spread out over 
a few degrees depending upon the initial phase of the orbit. 

The first of these results means that theoretical esti- 
mates, based on ir, of the expected fraction of accreted gas 
disk which will settle down to give polar rings will be too 
high. However the error is not likely to be greater than about 
10 percent, much lower than the uncertainty in the observed 
frequency of polar ring systems. 

More interesting is the second result which may pro- 
vide a mechanism for overcoming the angular momentum 
barrier which prevents accretion into the core of a galaxy. If 
material is accreted at an inclination between the measured 
maximum and minimum values of is, and if the accreted ma- 
terial is spread out over a range of phases, then precession 
will occur around both the minor and the major axes. Inter- 
actions between the two components can then lead to a large 
reduction in angular momentum. Indeed, once each compo- 
nent has precessed through 180 degrees then the planes of 
their orbits again coincide but their angular momenta are 



© 0000 RAS, MNRAS 000, 000-000 



4 Thomas, Vine & Pearce 



Figure 4. 
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